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A general analysis of Meissner effect and spin susceptibility of a uniform superconductor in an 
asymmetric two-component fermion system is presented in nonrelativistic field theory approach. 
We found that, the pairing mechanism dominates the magnetization property of superconductivity, 
and the asymmetry enhances the paramagnetism of the system. At the turning point from BCS to 
breached pairing superconductivity, the Meissner mass squared and spin susceptibility are divergent 
at zero temperature. In the breached pairing state induced by chemical potential difference and mass 
difference between the two kinds of fermions, the system goes from paramagnetism to diamagnetism, 
when the mass ratio of the two species increases. 
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It is well-known that there are two fundamental features of an electromagnetic superconductor, the zero resistance 
qq | and the perfect diamagnetism. The latter is also called Meissner effect The key quantity describing the Meissner 
effect is the Meissner mass or penetration depth. In the language of gauge field theory, the Meissner mass is the 
mass of the electromagnetic field obtained through the spontaneous breaking of local U(l) gauge symmetry, i.e., the 
Anderson-Higgs mechanism 2] . Recently, the study on superconductivity is extended to color SU (3) gauge field of 
Quantum Chromo dynamics (QCD) at finite temperature and baryon density Q. 

In the linear response theory, the Meissner effect is defined in the static and long wave limit u> — > 0, q — > of the 
■ external magnetic potential A(q). From the microscopic BCS theory the electric current density can be expressed as 

o: Q 
in ■ 

jq) = -Aq 1 + 7T^ / dpp 4 ^^) , 1 

where m,e,p and n are respectively the mass, electric charge, momentum and density of electrons, £a = 
yj (p 2 /2m — /j,) 2 + A 2 is the quasi-particle energy with electric chemical potential fj, and energy gap A, and f(x) 
the fermion distribution function. Since df(e<\)/de<\ < 0, the second term in the bracket on the right hand side is a 
• • . paramagnetic one and cancels partially the diamagnetism characterized by the first term. However, the total Meissner 
mass squared keeps positive in normal superconductor with BCS pairing mechanism. At zero temperature, due to the 
limit df(eA)/de& — > — S(ej\), the second term in disappears automatically, and there is no paramagnetic part. In 
addition to the perfect diamagnetism, the Meissner effect includes also the property of magnetic flux expulsion upon 
cooling through the critical temperature corresponding to the thermodynamic critical field. 

Another quantity to describe the magnetization property of a superconductor is the spin susceptibility. Since 
an electron carries a Bohr magneton, a cold free electron gas exhibits Pauli paramagnetism^]. However, at zero 
temperature the spin susceptibility x of a metallic superconductor is zero, it does not possess Pauli paramagnetism. 
The physical picture is clear: The two electrons in a Cooper pair carry opposite spin. From the microscopic BCS 
theory, the spin susceptibility of a superconductor can be written as[l| 

X - 2 £F r ^(- d 4^-) , (2) 



\p 37T 2 n J \ de A 

where xp is the Pauli susceptibility of a normal electron gas, and cf the electron energy at the Fermi surface. Due 
to the above mentioned limit of df(ej\)/dej\, the spin susceptibility of BCS type superconductor is zero at T = 0. At 
finite temperature, x is nonzero because of the thermo excitation in the superconductor. 

The above discussed diamagnetic Meissner mass and zero spin susceptibility are only for normal BCS superconductor 
where the two fermions participating in a Cooper pair are symmetric, i.e., they have the same chemical potential, 
the same mass, and in turn the same Fermi surface. In many physical cases, however, the difference in chemical 
potentials, or number densities, or masses of the two kinds of fermions results in mismatched Fermi surfaces. Such 
physical systems can be realized in, for instance, a superconductor in an external magnetic field£| or a strong spin- 
exchange field [EIEI3' an electronic gas with two species of electrons from different bands 8], a superconductor with 
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overlapping bands0,0|, a system of trapped ions with dipolar interactions a mixture of two species of fermionic 
cold atoms with different densities and/or masses \§L Il2j. an isospin asymmetric nuclear matter with proton-neutron 
pairing^, and a neutral quark matter in dense QCDjll M IlE E3, EE ISO, EL E3, E3, El, In the study 
on superconductivity in an external magnetic field, Sarma |j| found an interesting spatial uniform state where there 
exist gapless modes. However, compared with the fully gapped BCS state, the Sarma state is energetically unfavored 
and therefore instable. A spatial non-uniform ground state where the order parameter has crystalline structure was 
also proposed for such type of superconductors by Fulde and Ferrell and Larkin and Ovchinnikov, the so-called LOFF 
state. In this ground state, the translational and the rotational symmetries of the system are spontaneously broken. 
Recently, the above spatial uniform ground state prompted new interest due to the work of Liu and Wilczekjg. They 
considered a system of two species of fermions with a large mass difference. The stability of the state has been 
discussed in many papersj^ 0, HE HH, H3- It is now accepted that the Sarma instability can be avoided by two 
possible ways, finite difference in number densities of the two species |28[ or a proper momentum structure of the 
attractive interaction between fermions|23- In these states, the dispersion relation of one branch of the quasi-particles 
has two zero points at momenta p\ and p2, and at these two points it needs no energy for quasi-particle excitations. 
The superfluid Fermi liquid phase in the regions p < pi and p > p2 is breached by a normal Fermi liquid phase in the 
region p\ < p < p2- The temperature behavior of such a breached pairing (BP) state is very different from that of a 
BCS state, the temperature corresponding to the maximum gap is not zero but finite [la. 124, l29| . 

Since the dispersion relation controls the Meissner mass and the spin susceptibility, as shown above, it is natural 
to guess that the change in in breached pairing superconductor will modify the Meissner effect and spin sus- 
ceptibility significantly. Recently, it is found that the Meissner mass squared of some gluons in two flavor neutral 
color superconductor are negative [3(| On], which indicates that the quark matter in breached pairing state exhibits 
a paramagnetic Meissner effect (PME). In condensed matter physics, PME was observed first in high temperature 
superconductors such as ceramic samples of BioSroCaCuiOfi-LxpA. l33|. and later in conventional superconductors 
such as Nb|34L 1351 lipl l37| . It was suggested that the paramagnetic response might be a manifestation of d-wave 
superconductivity However, it seems that it is not necessarily to do anything with the d-wave analysis for the 
PME in conventional superconductors [35|. It is now widely accepted that the PME in these materials is most likely 
due to extrinsic mesoscopic or nanoscale disorder 39J. In this paper, we will investigate the Meissner effect and spin 
susceptibility in an asymmetric two-component fcrmion system in nonrelativistic case, and try to prove that the Meiss- 
ner effect and spin susceptibility of a superconductor is dominated by the pairing mechanism, and the paramagnetism 
and nonzero spin susceptibility are universal phenomena of superconductors with mismatched Fermi surfaces. We will 
model the pairing interaction by a four-fermion point coupling, which is appropriate for both electronic system, cold 
fermionic atom gas, nuclear matter and dense quark matter. Since our purpose is a general analysis for the Meissner 
effect and magnetization property, we will neglect the inner structures of fermions like spin, isospin, flavor, and color, 
which are important and bring much abundance while are not central for pairing. 

The paper is organized as follows. In Section [HJ we review the BCS theory in a symmetric fermion system and 
show how to calculate the Meissner mass and spin susceptibility in a nonrelativistic field theory approach. In section 
IIIII we extend the investigation to an asymmetric two-component fermion system with mismatched Fermi surfaces 
and derive the universal formula of Meissner mass squared and spin susceptibility. We then consider two kinds of 
mismatched Fermi surfaces induced by chemical potential difference and mass difference in Sections IIVI and In 
Section IVII we apply our general discussion to a relativistic system with spin structure and, as an example, reobtain 
the 8th gluon Meissner mass in neutral color superconductor. We summarize in Section IVIII We use the natural unit 
of c = h = 1 through the paper. 



II. SYMMETRIC FERMION SYSTEM 



In this section we review the BCS theory, the Meissner effect and the spin susceptibility in a symmetric fermion 
system in a field theory approach. We start with a system containing two species of fermions represented by a and b, 
described by the following nonrelativistic Lagrangian density with a four-fermion interaction, 

/ d V 2 \ 

C= }^ ipi(x) {-Q- + y~ + Mj V'iO) + gtp a (x)ipb(x)*l>b(x)*l>a(x) , (3) 

i—a,b ' 

where i/>(x), tp{x) are fermion fields for the two species, and the coupling constant g is positive to keep the interaction 
attractive. For the symmetric system, the two species have the same mass m and chemical potential fi. 
The key quantity to describe a thermodynamic system is the partition function which can be defined as 

Z= f[d^ a }[dMdMd^ dT ^ d3xC (4) 
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in the imaginary time (r) formulism of finite temperature field theory. According to the standard BCS approach, we 
introduce the order parameter A(x) of superconductivity phase transition and its complex conjugate A*(x), 

A(x) = g{^ b {x)^ a {x)) , A*(x) = g(Mx)Mx)) , (5) 

where the symbol ( ) means ensemble average. Since we focus in this paper on uniform and isotropic superconductor, 
we take the condensate to be ^-independent and real in the following. Introducing the Nambu-Gorkov space|lj defined 
as 

*=(t:) , *=(^« <M , (6) 



the partition function in mean field approximation can be written as 

Z MF = [[dnd*W< lT f d3 <*S- 1 *-W 2 / a ) (7) 



with the inverse of the mean field fermion propagator 



S ~'= ( A- " -£-£-„) ■ M 

\ or 2m ' / 



Taking the Gaussian integration in path integral Q and then the Fourier transformation, the thermodynamic potential 
of the system can be expressed as 

n = y - T E/(0 Trl ^"Vn,P) (9) 

in momentum space, where ^2 n is the fermion frequency summation in the imaginary time formulism. The first term 
is the mean field contribution, and the second term comes from the quasi-particle excitations with the inverse of the 
propagator 

s-Wl = f^ e % * ) do) 



in terms of momentum and frequency w„ = (2n + l)irT, where e p — — /i is the fermion energy. 

To determine the order parameter, the occupation number of fermions, the Meissner mass, and the spin susceptibility 
as functions of temperature and chemical potential, we need to know the fermion propagator itself. Using matrix 
technics it can be easily evaluated as 

rt- ~\ ( 0il(^njP) 0i2(«w n ,p)\ , , 

with the elements 

^ = iu„ + ( P = iun - e P g^ = -A g^ = -A 

{iuj n ) 2 - e A ' {iuj n ) 2 - e 2 A ' (iuJ n ) 2 — e A ' " (iw„) 2 - e 2 A 



where £a = y £p + A 2 is the quasi-particle energy. The excitation spectra u>± (p) can be read directly from the poles 
of the fermion propagator, 

w±(p*) = ±e A . (13) 

It is easy to see that these excitations are all gapped with minimal excitation energy A. 

The fermion occupation numbers n a = (ip£ip a ),rib = {^t^Ph) can be either calculated from the derivative of the 
thermodynamic potential with respect to the chemical potential, or equivalently, obtained directly from the diagonal 
elements of the fermion propagator matrix, 

n a (p) = T^g n (iu n ,p) , n b (p) = -T^2g 2 2(iuJ n ,p) . (14) 
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After the Matsubara frequency summation, one has 

n a (p) = n b (p) = ~ ( 1 - + %( £A ) . (15) 
At zero temperature, the fcrmion distribution function /(€a) goes to zero, the occupation numbers are reduced to 



n a (p) = n b (p) = - [I - . (16) 

The gap equation which determines the gap parameter A as a function of T and /i self-consistently can be expressed 
in terms of the non-diagonal elements of the fermion propagator matrix, 

A -^E/(03^(--P) ; (17) 
which is equivalent to the minimum of the thermodynamic potential, 

After the Matsubara frequency summation, the gap equation reads 

A(l - gI A ) = (19) 

with the function 

1 f d 3 p l-2/(e A ) 

It is easy to see that there are two solutions of the gap equation i|19|) : One is A = which describes the symmetry 
phase, and the other is A ^ determined by 1 — gl& — which characterizes the symmetry breaking phase. 



A. Meissner Effect 

We show now how to calculate the Meissner mass in terms of the thermodynamic potential. Suppose the fermion 
field carries electric charge e and couples to a magnetic potential A. In mean field approximation, the magnetic 
potential is treated as an external and static potential, and the thermodynamic potential Q(A) of the system can be 
expanded in powers of A, 

fl(A) - n(0) + ^MfjAiAj + . . . , (21) 



with the coefficient 



M , = 9HKA) 
v dAidAj 



, ■ (22) 

A=0 



Since the thermodynamic potential is just the effective potential of the field system, the coefficients Mh can be defined 
as the components of the Meissner mass squared tensor. If the ground state of the system is isotropic, one has My = 
for i =/= j and M\ x = M| 2 = M| 3 , and the Meissner mass squared M 2 can be defined asplfl Elj 

(23) 



i o 2 n(A) 



M 



3 f-' dAidAi 



A 



In our model of four- fermion point interaction (J3J, the thermodynamic potential in the presence of external and 
static magnetic potential A can be expressed as 

A 2 f d 3 p 



"(A) = — -TY j ) J^T^A^^P) (24) 
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where the A-dependent propagator is defined as 

*-( iu; > 

A iuj n + e 

with the fcrmion energies = ^^m^ ^° cx t rac t the Meissner mass squared, we expand the propagator G^ 1 , 



g^ n ,p) = (^\ e p . * ) , ,25, 



Ga=S- 1 --P-A-^£-t 3 , (26) 
m 2m 

and its contribution to the thermodynamic potential, 

2 2 

Tr In g- 1 = TrlnQ- 1 - -p • ATrG - ^—A 2 Tr (0r 3 ) - -^(p • A) 2 Tr (gg) + ■■■ (27) 

m 2m 2mr 

in powers of A, where g(iu> n ,p) is the propagator matrix (|llf> in the absence of magnetic field. After the momentum 
integration, the linear term in A vanishes, and the Meissner mass squared M 2 can be read from the coefficient of the 
quadratic term in A 2 of the thermodynamic potential fl(A), 

M 2 = — E / 7^4 (&i Q22) + / 7?4t + ^ 22 + 2 ^ 21 ) ■ ( 28 ) 

m J (27r) J m z ,/ (2TT) ,i 3 

From the comparison with the fermion occupation numbers (|14|l . the first term on the right hand side is proportional 
to the total number density n, 

m ' / 27r) ci m 

with the definition 

^—3 (n a (p) + n 6 (p)) . (30) 

Employing the Matsubara frequency summation ^2 n GG calculated in Appendix lAl for the second term, we recover 
the well-known Meissner mass squared shown in text books 0, 

^.^(,+ 1^.1^ (31 ) 



m \ 3n 2 mn J 9ea / m 

with n, defined as 



1 r„,f«. (32) 



37r 2 m Jo V ^ e A 

The effective density 7i s is positive at any temperature and chemical potential, which means diamagnetic supercon- 
ductor for any symmetric fermion system. At low temperature limit and in the approach to the phase transition line 
of superconductivity, n s behaviors aspj 



*± = i-./^e-Ao/^ T ^ Q 
n V T 

n T 

- = 2(1--), T^T C , (33) 
n T c 

where Ao is the order parameter calculated by the gap equation 119fl at zero temperature, and T c the critical tem- 
perature determined by 1 — gIo(T c ) — 0. It is also necessary to note that the Meissner mass squared l|31|l satisfies the 
renormalization condition 

M 2 (A = 0) = 0. (34) 



6 



B. Spin Susceptibility 

Assuming the thermodynamic potential in the presence of a constant magnetic field B to be f2(-B), the magnetic 
moment M. and the spin susceptibility \ °f the system are defined as 



M 



dQ(B) _ d 2 n(B) 



dB 



b=o ' X dB 2 



(35) 

B=0 



In our model, the thermodynamic potential fl(B) in mean field approximation can be expressed as 

V(B) = ^--TJ2 J (0 Trlll ^ 1 (^n,p) , (36) 
where the independent propagator is defined as 

with the spin-up and spin-down fermion energies ej, = — fi + HoB and = ^ — /x — ^to-B, where /j,o is some 
elementary magneton such as the Bohr magneton or the nucleon magneton. 

To extract the magnetic moment and spin susceptibility from the expansion of f2(_B) in powers of B, we take the 
similar way used for the discussion of Meissner effect in the last subsection. We expand the propagator 

Qa 1 = S- 1 -fMiB , (38) 

and its contribution to the thermodynamic potential 

Tr\nG B l =Tr\ng- 1 - ^BTrQ -^{^B) 2 Tr{gg) + ■ ■ ■ (39) 

in powers of B. Substituting them into the thermodynamic potential, we obtain from the linear term in B the 
magnetic moment M. 



M 



(27T) 

and from the quadratic term the spin susceptibility^} 



^ T J2 J (&i + £ 22 ) = Vo(n b - n„) = , (40) 



= -m^E / ^(5nGn + 522 g22 + 2g 12 g 21 ) = -^lJ™dpp 2df{eA) 



Oca 



(41) 



where we have used again the Matsubara frequency summation of ^2 n GG shown in Appendix lAl It is easy to see 
that x = at T = and x = Xp = at T = T c . 

III. ASYMMETRIC FERMION SYSTEM 

We discuss now the fermion pairing mechanism, and the Meissner effect and spin susceptibility in an asymmetric 
fermion system, using the same approach for the symmetric system in Section [H] . Our asymmetric two-component 
system with different masses m a , mj, and different chemical potentials fi a , l^b is defined through the Lagrangian density 



i— a,fe 



a v 2 

t; l~ 7; 1- Mi V^W + 9i>a(xWb(z)i>b(z)'ipa(xi) • (42) 

ot zm,- 



The thermodynamic potential of the system in mean field approximation is just the same as © for the symmetric 
system, but the matrix elements of the fermion propagator in Nambu-Gorkov space are different, 

iuj n - e A + e s n iu! n ~ e A - e s 



t'll ~~r. To <T~ i y22 ~~r. y 

{lUn - eA) A - e A (lU n -eAr 



Gl2 = y. 75 o" ' ^21 = y. T5 g~ ' (^) 
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where eg and are denned as es = p 2 p , 6a = p 2 P with the fermion energies = ^ Ma and = — 



and ca — ye| + A 2 is t ne quasi-particle energy. The dispersion relations u>± (p) can be read from the poles of the 
fermion propagator, 

W±(p) = £A±£A ■ (44) 

Without losing generality we can choose > in the following. Different from the BCS mechanism for the symmetric 
fermion system, while one branch of the excitations lu + in the asymmetric system is always gapped, the other one w_ 
can cross the momentum axis and become gapless at the momenta p% and p 2 , where p\ and p 2 satisfy the equation 
w_ (p) = and p F < Pi < P2 < Pp with p F and p b F the Fermi momenta of the two species. 

The phenomena of gapless excitation is directly related to the breached pairing mechanism. The occupation numbers 
for fermions a and b defined in l|14f) become now 

MP) = \ (l - /(«_) + \ + 

n b (P) = ~\ (l + ^) /(*-) - ~ (l - ^) f(u>+) + 1 • (45) 
At zero temperature, they are reduced to 

n a (p) = ^1 - j^j 8(ca - e A ) , 

n b (p) = 1 - ~ + 0(e A - e A ) • (46) 

If there is no breached pairing, namely ea > £a, the two species have the same occupation number 

n a (p) = n b (p) = i (l - , (47) 

which comes back to the result (|16|l for m a = n%b and /i a = /i . However, when the breached pairing happens, the 
result 147|) is valid only in the momentum regions p < Pi and p > p 2 , and in the breached pairing pi < p < p 2 , we 
have 

n„(p) = , n 6 (p) = 1 . (48) 

In this case, the pairing between fermions is breached by the region p\ < p < p 2 , the system is in normal Fermi liquid 
state in the region p\ < p < p 2 and superconductivity state in the regions p < p\ and p > p 2 . 

For the asymmetric system, the gap parameter A is still determined through the self-consistent equation (|19f) . but 
the function 2a is changed to 



A. Meissner Effect 



Suppose the fermions a and b carry electric charges eQ a and eQb respectively, the thermodynamic potential in the 
presence of a magnetic potential A in mean field approximation is still in the form of l|24|l , but the propagator matrix 
is now a little bit different, 

g- A Hi^p)=(^S i0Jn + 4 -)=S- 1 -eT l p-A-^ (50) 
with the fermion energies = ^m' A ^ Ma, = ^Im^ ^ 6 > ana - ma trices Ti and T 2 defined as 

r ' = (H)' r *"(f (51) 



III b 
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Taking again the expansion of Tr In Q in powers of A, 



22 



Tr ln^ 1 = TrlnQ- 1 - ep • ATr (QT\) - — A 2 Tr{QY 2 ) - -(p ■ A) 2 Tr (SIWx) 

the Meissner mass squared M 2 can be extracted from the quadratic term in A of fi(A), 

M 2 = M 2 D + Mp , 
M 2 , T v- f d 3 p (Q 2 Ql 

The diamagnetic part is related to the number densities, 

M 2 D ^(^Ql + ^Ql\ e2; 

V m a m b J 

and the paramagnetic term can, with the help of the frequency summations in Appendix^ be expressed as 



(52) 



2QaQl 

m a m b 



Q\iQ: 



21 



(53) 



(54) 



M 2 = e 2 



d 3 p p 2 

lp = c 1 (2^y 

with the definitions 



Qa 



2 2 

u p v p 



f(w+) - f(0J-) 



Qa 2 i Qb 2 \ .', 

— U P + — V P f ( w ^ 
m a y m b t ' 



Qa 2 i Qb 2 \ .', , 

v + — u f (w_) 

m a F m b 1 ■ 



(55) 



■4-5 1 + 5 



2 / -i £ S 



= l-< 



(56) 



The first term in the square bracket of l|55() is a new term resulted fully from the asymmetric property of the 
system. Since f(u>+) < at any momentum, it is always negative. We see that the asymmetry between the 

paired fermions enhances the paramagnetism of the system. The second and third terms are negative due to the 
property f'(x) < for any x, they together will be reduced to the paramagnetic part of M 2 if we come back to 
the symmetric system. At zero temperature, we have the limit f(x) — > 9(— x) and f'(x) — » —8(x), the second term 
vanishes due to —d(uj+) — > 0, and the third term disappears only in normal superconductor but keeps negative in 
breached pairing state because of the property /'(w-) — > — S(ca — ea) 0- 

In the following two sections we will discuss in detail the paramagnetism of breached pairing superconductors 
induced by chemical potential difference and mass difference between the two paired fermions. Here we just point out 
a singularity of the Meissner mass squared at zero temperature in general case. At the turning point from gapped 
excitation to gapless excitation where the two roots p± and p 2 of ui-(p) = coincide, p\ — P2, the momentum 

integration of the third term of (|55|l goes to infinity, since / dpS(oj-) — > J dtu^ I ^( a - 1 -) ~~ * 00 due to 

at p± — P2, and then the total Meissner mass squared becomes negative infinity at this point. 

It is easy to check that there is still the renormalization condition for the total Meissner mass squared, 



dp 







M 2 (A = 0) = M 2 (T = T c ) = 



(57) 



at the critical temperature T c 



B. Spin Susceptibility 

As shown in Section [HI the magnetic moment M. and spin susceptibility x are, respectively, the coefficients of 
the linear and quadratic terms in the magnetic field B in the thermodynamic potential £!(£?). For the asymmetric 

system, the fermion energies in the independent propagator matrix Ij37|l are defined as = ^ A*a + HodaB, and 

2 

e p = 2m — ^o9bB, where g a and g b are constants related to the quantum numbers of angular momentum of 

species a and b. 

Expanding the propagator 

Gb 1 =Q- x -MoBr 3 , (58) 
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and 



Tr In Q B X = TrlnQ- 1 - ti BTr(GT 3 ) - ^ B) 2 Tr{gT 3 gT 3 ) + ■■■ (59) 



in powers of B with the matrix 



g b 

we obtain from the expansion of £l(B) the magnetic moment 

(27T) 



r 3 = ( 9 n ° ) > (60) 



= -HoT^2 / 7^3 (So^ll + 9bQ2l) = Mo(ff6«b - 5a™a) , (61) 

and the spin susceptibility 

(62) 



d 3 p 



-Mo 



/ \2 2 2 /( Ll '+) , / 2 1 2\ 2 /»/ / \ 1 / 2 1 2\ 2 j-// \ 

{9a -9b) u v ^{9aU +g b v ) f + (g a v 1 + g b u ) f (w_) 



(2tt) 3 

where we have again taken into account the frequency summations listed in Appendix ^ Similar to the discussion 
for the Meissner mass, the first term here is new and fully due to the asymmetry property between the two species, 
and x is divergent at zero temperature at the turning point from gapped to gapless excitations due to the behavior 
of the term with f'(ui-) in the limit T — ► 0. 

We now turn to the details of the Meissner effect and spin susceptibility in breached pairing state induced by 
chemical potential difference and mass difference between the two paired fermions. 

IV. ONLY CHEMICAL POTENTIAL DIFFERENCE 

We first discuss the Meissner effect and spin susceptibility induced by chemical potential difference only. It is 
convenient to replace the chemical potentials of the two species [i a and ^ by their average /2 and difference Sfi defined 
as 

[i-a + c Lib - Ma ,„„s 

M = 2 ' M = 2 ' ' ' 

Without losing generality, we can set Sfi > 0. With p, and (5/i, the dispersion relation of the elementary excitations 
can be written as 



w±(^)=*M±y(^-A) +A 2 . (64) 
It is easy to see that only under the constraint 

A < Sfj. , (65) 
there is breached pairing in the momentum region p\ < p < P2 with 



Pi = J 2m (/i - \J 8y? - A^J , p 2 = ^2m(^i+ V^A 



/i 2 - A 2 . (66) 



A. Meissner Effect 

For simplicity, we set Q a = Qb = 1. In this case the diamagnetic and paramagnetic parts of the Meissner mass 
squared take the form 



rie 



2 



Mi = — , (67) 
m 

2 poo 4 

M 2 P = — / dpf, [/'(e A - &H) + f(e A + 5»)] . 
m z J bit 1 
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At zero temperature, the paramagnetic part is evaluated as 



M 2 P 



e 2 [°° p 4 
— 2 / d P~^i S ( £A ~ ^) 



(68) 



It is easy to see that only in the breached pairing state, namely, A < Sfi, Mj, is nonzero. After the momentum 
integration, the total Meissner mass squared can be expressed as 



M = 



ii< - j , 5pi6(5n — A) 




with the parameter r/ defined as 



and the total fermion density 



Pl+pl 



A 2 



67r 2 n 



P1+P2 
6tt 2 



dpp 2 u 



2„.2 
P 



dpp 2 vl 



P-2 



(69) 



(70) 



(71) 



From the definition of u 2 and u 2 and their relation to the occupation numbers n a (p) and n&(p), the first and second 
integrations in (|71|l are, respectively, the upper and lower shadow regions of n a {p) and rty,(p) in three dimensional 
case, shown diagrammatically in Fig^ Since n a (pi) = nb(pi) > 1/2 and n a (j>2) — n\ } {p2) < 1/2, the contribution of 
the shadow regions to the total fermion density n is much smaller compared with the term (p\ +p^)/Q'K 2 and we have 
approximately, 



77 ~ 1 , M 2 ~ — I 1 - 

m 



SfiOiSfi - A) 
TP 



A 2 



(72) 



which means global paramagnetism in the asymmetric fermion system, if the breached pairing happens. 




P| Pi P 2 PP 

FIG. 1: The schematic occupation numbers n a and ri(, as functions of momentum, an d Pf are Fermi momenta of species 
a and b, and pi and P2 are the two roots of the dispersion equation w_(p) = 0. The pairing state is breached by the normal 
fermion liquid state in the region pi < p < P2- 



B. Spin Susceptibility 



We take g a = gb = 1 for simplicity. While the magnetic moment M. disappears automatically for normal supercon- 
ductor, it is no longer zero in the breached pairing superconductor, 



M = Ho(n b - n a ) = /_i 



pl 



p\ 



6tt 2 



(73) 



11 



The physical picture is clear: while the paired fermions in the region p < p% and p > pi have no contribution to the 
magnetic moment, the unpaired fermions in the region pi < p < p 2 do contribute to M. 
The spin susceptibility in this case is reduced from to 



X = -Mo 



dp^ (f'(e A - Sfi) + f'(e A + Sfi)) . (74) 



At zero temperature, it is nonzero only in the breached pairing state characterized by A < S/i, 

X = 2 R „ = A Pi +P2) ■ (75) 
2tt 2 y'Sp, 2 - A 2 

From the comparison with the Pauli susceptibility xp, we have the relation 

X Sfi 



Xp ^5 



(76) 



V. BOTH CHEMICAL POTENTIAL DIFFERENCE AND MASS DIFFERENCE 

We now consider the magnetization property of the superconductor with both chemical potential difference and 
mass difference between the paired fermions. To simplify the calculation, we still take Q a = Qb = 1 and g a = Qb = 1- 
From the dispersion relations 



w±(p) 



V 



■2,,,, ■ WG^ - '' 1 +A2 



(77) 



with the reduced masses tua = 2T " a "' b and ms = 2m l" b , the condition for the system to be in breached pairing 
state is 



A < A, = 



|m b /Ub - m a u a \ _ \Xfi b - Mai 



2y/m a m b 



2y/X 



(78) 



where A = mb/m a is the mass ratio, and the corresponding region of breached pairing is located at p\ < p < pi in 
momentum space with 



Pi = 



\Jm a {p, a + A^b) - V (pa - Vb) 2 - 4AA 2 



P2 = \ m a {fl a + A^if,) + y/ (p a - XfJ, b ) 2 - 4AA 2 



(79) 



A. Meissner Effect 



We now calculate analytically the Meissner mass squared at zero temperature. When there is no breached pairing, 
we have from the general expressions (|53|l to l|55|) . 



( A + l)ne 2 _ (A-l) 2 e 2 f°° ^ p 4 u 2 p v 2 p 
2Xm a X 2 m 2 J 6ir 2 ca 



.'. - V \ 2 % dp^^JL. (80) 



In the state with breached pairing, taking into account the relation 

d -^-p(^-^)-^(ul~H) (") 

op \ nib ma I Xm a F F/ 

for the integrated function with /'(oj_) in (|55|l . the total Meissner mass squared in the breached pairing superconductor 
can be written as 
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M 



f2 



6Am a 7r 2 
3(A + 1) 



Pi 



[i + (a-iK] j 

|l-(A + l)«f| 



Pi 



[1 + (A-1)^ 2 2 ] 2 
|1-(A + 1>!| 



in 



dpp 2 ul 



l>2 



dpp 2 v 2 - 2(A - 1) 




(82) 



with the shorthand notations v\ 



2 

Pi allU "2 - -p 2 - 

It is easy to check that for A = 1, we recover the result obtained in Scction lTVl for the case with only chemical potential 
difference. For A > 1, we again take the approximation of neglecting the integration terms in (|82p. Considering the 
relation 



(i + (x-i)v 2 2 y 



\l-{\ + l)v 2 \ 



for any A, the sign of M 2 depends on the quantity 



f}{\,v{) = \- 



(1 + (A-1W) 2 
|1-(A + 1K| 



(83) 



(84) 



Fig|21shows (3 as a function of v\ at A = 2 and A = 10. Taking into account the condition v\ = n a (pi) — rn,(pi) > 1/2, 
we focus on the behavior of (3 in the region of v 2 > 1/2. At A = 2, f3 is negative in this region and results in 
negative Meissner mass squared and paramagnetism of the system, which is continued with the conclusion in Section 
IIVI However, at A = 10, (3 becomes positive in the interesting region, the Meissner mass squared tends to be 
positive and the system tends to be diamagnetism. In fact, for very large A, v\ becomes extremely small 12J, and 

1 — (l + (A — l)v|) /|1 — (A + l)w|| ~ 0, the system contains approximately the diamagnetic term only. 




FIG. 2: The parameter (3 as a function of Vi . The dashed and solid lines correspond to A = 2 and A = 10, respectively. 



B. Spin Susceptibility 

In the case with both chemical potential difference and mass difference, the magnetic moment M. in breached 
pairing state takes still the form of (|73l) , but pi and P2 determined by the dispersion relation cj_ (p) = are shown in 
(|7§|l . As for the spin susceptibility x, its general expression at finite temperature is still l|74|). At zero temperature, 
it is reduced to 



Am a Mo ( Pi P2 



2tt 2 \\u 2 ~Xv 2 \ \u 2 -\v 2 \ 



(85) 
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VI. EXTENSION TO RELATIVISTIC SYSTEMS 



We have investigated the Meissner effect and spin susceptibility in an asymmetric system of two kinds of fcrmions 
with different chemical potentials, masses, charges, and magnetic moments in nonrelativistic case. We found that the 
magnetization property of breached pairing superconductor is very different from that of BCS superconductor, and 
the system tends to be more paramagnetic. What is the situation in relativistic case? Are these exotic phenomena 
just a consequence of nonrelativistic kinematics? In this section, we extend our discussion to relativistic systems, 
which is relevant for the study of color superconductivity at high baryon densitvp^. E^. EEL Effll l3ll| . We will see 
that there is still paramagnetic Meissner effect in breached pairing state in relativistic superconductors. 



A. Without Dirac Structure 



As a naive calculation, we first neglect the antiparticles and take the same Lagrangian density (|42|l . The relativistic 
effect is only reflected in the fermion energies e° = \/p 2 + m 2 — /.i a , e b p = \/p 2 + m 2 — /.tj,. The formulas for Meissner 

mass squared and spin susceptibility in Section IlIII still hold if we replace rrnii = a, b) there by \/p 2 + mf. As an 
example, we list here the diamagnetic and paramagnetic parts of the Meissner mass squared in the case with only 
chemical potential difference between the two species, 

Ml = ^ , (86) 
m 

e 2 f°° d 4 
M 2 P = — / d V ^—^ [/' (e A - 5n) + f (e A + Sp)} 

QTT Z Jq p Z + Tfl 

with the relativistic fermion density 

-J 7 =Kp )+n b (p)} . 87 
( 27r ) vp + m 

In ultra relativistic limit m — > and at zero temperature, they can be evaluated as 

eV A Sfi 2 - A 2 



M d - ^(1 + ^—^) , 



1 " —\^^^) ' 

P 3tt 2 ^/Sj? - A 2 V M 

where in the calculation of we have taken again the approximation used in deriving (|72|l . The paramagnetic part 
is automatically zero in normal superconductor with A > Sfi, but negative in breached pairing superconductor with 
A < S/i. Putting the two terms together, the total Meissner mass squared can be expressed as 

e 2 /Z 2 / 5y 2 -A 2 \ ( Sj^Sp-A)" 
3tt 2 ^ + fl 2 ){ vV - A 2 



^-^1+^ 3- 7 -M . (89) 



It is negative in the case of A < ^8/96^ < 8fi. Therefore, a relativistic breached pairing superconductor is also 
paramagnetic. 



B. With Dirac Structure 



We study now a more realistic relativistic model containing two kinds of fcrmions. The Lagrangian of the system 
is defined as 

a—a,b a.J3—a,b a.J3—a.b 

where 7 M = (7 , 7 1 , 7 2 , 7 3 ) = (/3, (3a 1 , /3a 2 , /3a 3 ) and 7 s = i7°7 1 7 2 7 3 are Dirac matrices with f3 and a 1 being an- 
ticommuting matrices, /3 2 = a 2 = l,{oti,otj} — for i ^ j, and {ai,/3} = 0, their covariant counterparts 7^ and 
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75 are defined as 7 M = (7 , — 7 1 , — 7 2 , — 7 3 ) and 75 = 7 s , V and tp are Dirac spinors, ip c = Cip T and t/> c = ip T C 
are charge-conjugate spinors, C = i~f 2 j° is the charge conjugation matrix, the superscript T denotes transposition 
operation, and T\ is the first Pauli matrix with the elements Tf a = r\ h — and rf 6 = t\°- = 1 . 
For convenience we define the Nambu-Gorkov spinors 



, * = ( i>b ^6 ) • 



(91) 



Introducing the order parameter 



a : f3—a 7 b 



(92) 



and taking it to be real, the thermodynamic potential of the system in mean field approximation is still in the form 
of © with the propagator matrix defined in the 4-dimensional Nambu-Gorkov space, 



g- 1 



z 75 A \ 

i Tj a [g^ 1 

z 75 A 

V ?75 A [So"]- 1 J 



(93) 



where So is the free propagator, 



[Go 



±1-1 



«w„7o - 7 • p - m a ± HaJo 



(94) 



with the vector matrix 7 = (7 , 7 2 , 7 3 ). 

Assuming that the fermion field ip a with charge e can couple to some U(l) gauge potential A^, the Meissner mass 
squared can be extracted from the expansion of the thermodynamic potential in powers of the external magnetic 
potential A. By separating the A-dependent propagator (|93|l into two terms, 



with the self-energy matrix E defined by 



g a (iuj n ,p) = g 1 + ■ A 



10 
0-100 
10 
-1 



(95) 



(96) 



and taking the coefficient of the quadratic term in A of f2(A), the Meissner mass squared can be evaluated as 

r d 3 p 3 r •" 
M 2 = e 2 T^2 I f 2 P" ^ Tr +0227*0227* +0337*^337* + GiAl l G4Al % - 20127*0217* - 20347*0437* • (97) 

In the case of m a — rrib, it is straightforward to write down explicitly the propagator matrix elements Sij, by 
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employing the method used in |42t |43| , 

GzA{iu ni p) 
Gi3{iu ni p) 
with the quasi-particle energies 

and the energy projectors 



i0Jn 


+ c p 


— 5fi 


(iuj n - 




- fe) 2 


lUJn 


- e p 


— 


(lU n - 


5fx) 2 


- (e A ) 2 




iA 




(iuj n - 




- fe) 2 




iA 




{iu n - 


Sp) 2 


- fe) 2 


lUJn 


+ e- 


+ Sfi 


(iu n + 




- (e A ) 2 


iuj n 


~ 6 p 


+ Sn 


(iuj n + 


sri 2 


(e A ) 2 




iA 




(iLO n + 


Sri 2 


- (e A ) 2 




iA 




(icj n + 


Sri 2 


- (e A ) 2 



A+7o + 
A_7o + 
A-75 + 
A+75 + 
A+7o + 
A-70 + 77 



A+75 





-4- 


5[i 


(iuj n - 
iuj n 


6ri 2 - 

+ 4- 


(4) 2 


(ioj n - 


sri 2 - 

iA 


(4) 2 


(iuj n - 


sri 2 - 

iA 


(4) 2 


(iu n - 


sri 2 - 

-4 + 


(4) 2 

S/jL 


(iu> n + 
iuj n 


sri 2 - 
+ 4 + 


(4) 2 


(iuj n + 


sri 2 - 

iA 


(4) 2 


(iuj n + 


Sri 2 ~ 
iA 


(4) 2 


(iU n + 


sri 2 - 


(4) a 


J(4) 2 


+ A 2 , 





A-70 
A+7o 
A+75 
A_75 



A+7o 
A+75 



A± 



- ( 1 ± 7o ^ ' p 



Vp 2 + 



(98) 



(99) 



(100) 



It is easy to see that in the ultra relativistic limit m — > 0, the expression l|97|l is just the same as the 8th gluon's 
Meissner mass squared in the two flavor gapless color superconductor [3(], El, namely, 



M 



2e 2 fl 2 
3tt 2 



6fi6(5(i - A) 



(101) 



Since we did not consider here the non-Abelian structure of the color superconductivity, the negative Meissner mass 
squared for the 8th gluon is just a reflection of the breached pairing mechanism. 



VII. SUMMARY 



We have investigated the relation between the pairing mechanism and magnetization property of superconductivity 
in an asymmetric two-component fermion system coupled to a magnetic potential. In the frame of field theory 
approach, we derived the dependence of the Meissner mass squared, magnetic moment, and spin susceptibility of 
the system on the chemical potential difference, mass difference, charge difference, and magnetic moment difference 
between the two kinds of fermions. Compared with the superconductor formed in a symmetric system where the 
Meissner mass squared is globally diamagnetic, there is no magnetic moment, and the spin susceptibility disappears 
at zero temperature, we found the following new magnetization properties for the asymmetric system with mismatched 
Fermi surfaces between the paired fermions: 

1) The asymmetry leads to a new paramagnetic term in the Meissner mass squared and a new term in the spin 
susceptibility, and the magnetic moment is no longer zero in asymmetric systems. Note that the new terms and 
the finite magnetic moment do not depend on the pairing mechanism, they are only the consequence of asymmetry 
between the two species. 

2) At the turning point from BCS to breached pairing state, the Meissner mass squared and spin susceptibility are 
divergent at zero temperature. 
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3) In the breached pairing state induced by chemical potential difference between the two paired fermions, the 
Meissner mass squared is paramagnetic at zero temperature. 

4) In the breached pairing state with not only chemical potential difference but also mass difference between the two 
kinds of fermions, the system at zero temperature is paramagnetic at small mass ratio and tends to be diamagnetic 
when the ratio is large enough. 

While the paramagnetic Meissner effect and finite spin susceptibility discussed above are interesting, how to under- 
stand them correctl y an d their reflection on physically observable quantities are not clear. By comparin g th e BP and 
LOFF states|44j,|45|,|4g, the paramagnetic Meissner effect might be a signal of instability of the BP state|44j,|45j. The 
thermodynamic potential f2(Q) of a LOFF state with a nonzero momentum 2Q of the Cooper pair can be obtained 
by replacing the magnetic potential eA in the thermodynamic potential f2(A) derived above by the momentum Q. 
If the uniform BP state is the ground state, the thermodynamic potential must be the minimum at Q = 0, namely, 
dfl/dQi = and k = d 2 ft/dQ 2 > at Q = 0. Through the obvious relation M 2 = e 2 n, negative Meissner mass 
squared leads automatically to negative n. Therefore, the paramagnetic Meissner effect may be a signal that the 
LOFF state is more favored than the BP state. However, the above argument is obtained from the study for systems 
with fixed chemical potentials, it is not clear if it is still true for systems with fixed number densities of the two 
species. As we know, the existence of the LOFF phase in conventional superconductors has still not been convincingly 
demonstrated in an y m aterial. If the above argument is true, the LOFF state might be observed in trapped atomic 
fermion systems [471 l48j |. Our research in this direction is in progress. 
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APPENDIX A: FERMION FREQUENCY SUMMATIONS 



We calculate in this Appendix the fermion frequency summations TJ2 n GnGn, TJ2 n G22G22, an d G12G21 
in the Meissner mass squared and spin susceptibility in Sections IIIII IIVI and [V] From the decomposition of these 
summations, 

r5>x0 u = tJ2 rn~~ '"L = A + (^A + e s ) 2 B + 2e s C , 



[(iu> n - 




-el? 


(iu> n - 


- e A - 




[(iu n - 


tA) 2 


~e\} 2 




A 2 





[{iuj n - eA) 2 - e 2 A } 2 



A+(e A ~e s ) 2 B~2e s C , 



A B 



with the definitions 



A = TV — ^ 2-, B = Ty— ; 

^ (lUJ n - 6 A ) 2 - 6 2 A - , , 



[(iu) n - e A ) 2 - 



2 _ c 2 12 



-AJ 



{iu) n - e A + £A) 2 {iu n - e A - ea) 
we need to complete the summations A, B and C only, 



A = 
C = 



/(-a;,) + /(«;+)- 1 ^ B = 

2eA 2eA de& 

f(-LJ-)+f(u+)-l 1 df(-u-) 

44 



2e A 



2eA deA 

Defining the function f'(x) = df(x)/dx, we can express the summations we need as 

T^GnGn = u p v p \- v p f (cj_) + Upf 

T 2^ G22G22 = u p v p + U p f + Vpf (cj+) , 



£A 



^^^12^21 



-u p v p 
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where the energies £a,^s,^a, the dispersion relations ui+,ui-, and the functions Up,Vp are defined in Section IlIII 
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